Forecasting Sauter mean droplet size and examining the range of droplet sizes in a Tenova liquid–liquid extraction column

A new type of Tenova pulsed extraction column was introduced in 2017. It is the newest generation of pulsed columns. Due to the internal equipment of this column and the lack of moving parts and the simplicity and speed of repairs and maintenance, it has been the focus of researchers in recent years. No correlations for predicting the mean drop size and drop size distribution of the Tenova column have been reported. The Sauter mean drop diameter and drop size distribution are investigated for a Tenova pulsed column with a diameter and an active height of 7.4 and 73 cm, respectively. Three standard chemical systems of isobutyl acetate-water, isobutanol-water, and toluene-water have been used. The effects of pulse intensity, dispersed and continuous phase flow rates have been taken into account. In each experiment, 200–300 drops have been analyzed in a total of 10,000 drops. The investigation covered a spectrum of physical properties, notably surface tension (within a range of 1.75–36 mN/m). Operating conditions including pulse intensity (in the range of 0.2–2 cm/s) and the flow rate of continuous and dispersed phases (in the range of 8–30 L/h) have been investigated. Methods based on artificial intelligence (AI) such as multilayer perceptron neural networks and gene expression programming were combined with a dimensional analysis approach to provide a new approach to estimating the mean drop diameter (d32). Experimental results have been compared with the equations found by other researchers in similar columns. The variation of drop size distribution has also been experimentally obtained.Methods based on artificial intelligence (AI) such as multilayer perceptron neural networks and gene expression programming were combined with a dimensional analysis approach to provide a new approach to estimating the mean drop diameter (d32). Experimental results have been compared with the equations found by other researchers in similar columns. The variation of drop size distribution has also been experimentally obtained.

No correlations for predicting the mean drop size and drop size distribution of the Tenova column have been reported.The Sauter mean drop diameter and drop size distribution are investigated for a Tenova pulsed column with a diameter and an active height of 7.4 and 73 cm, respectively.Three standard chemical systems of isobutyl acetate-water, isobutanol-water, and toluene-water have been used.The effects of pulse intensity, dispersed and continuous phase flow rates have been taken into account.In each experiment, 200-300 drops have been analyzed in a total of 10,000 drops.The investigation covered a spectrum of physical properties, notably surface tension (within a range of 1.75-36 mN/m).Operating conditions including pulse intensity (in the range of 0.2-2 cm/s) and the flow rate of continuous and dispersed phases (in the range of 8-30 L/h) have been investigated.Methods based on artificial intelligence (AI) such as multilayer perceptron neural networks and gene expression programming were combined with a dimensional analysis approach to provide a new approach to estimating the mean drop diameter (d 32 ).Experimental results have been compared with the equations found by other researchers in similar columns.The variation of drop size distribution has also been experimentally obtained.Methods based on artificial intelligence (AI) such as multilayer perceptron neural networks and gene expression programming were combined with a dimensional analysis approach to provide a new approach to estimating the mean drop diameter (d 32 ).Experimental results have been compared with the equations found by other researchers in similar columns.The variation of drop size distribution has also been experimentally obtained.Liquid-liquid extraction is one of the most widely used separation processes 1,2 .This method is the most common separation process in nuclear and hydrometallurgy industries and the production of relatively pure chemical compounds and environmental waste treatment [3][4][5][6][7][8] .Pulsed extraction columns are efficient liquid-liquid contactors that have been designed and studied in different types of internal equipment such as packed, tray, and disc and doughnut.One of the advantages of pulsed columns is the absence of moving parts and simplicity and the speed of repairs and maintenance of such columns, while the high efficiency of such columns has caused special attention to such columns [9][10][11][12][13][14] .The new generation of disc and doughnut columns has been introduced as "Tenova column".Tenova's internal equipment has new contact factors on disc and doughnut plate alignments and optimized spacing between plates to achieve less reverse mixing, higher column dispersed phase inventory, and improved mass transfer with a uniform flow rate during operation.The design of an extraction column requires the examination of the separation performance and the detailed examination of the hydrodynamic parameters such as the drop size and the amount of the dispersed phase as well as the flooding velocities to determine the column capacity.

List of symbols
The aim of this study was to investigate the hydrodynamic behavior of Tenva pulse columns and achieve changes in the Sauter mean drop diameter and drop size distribution along the column as the basis for the design of operational units.Experiments were performed along the column using three different liquid-liquid systems.
We incorporated this information to highlight the novelty of the Tenova column and the absence of predictive relationships for crucial parameters, specifically, the average diameter and distribution of drops essential in designing these columns.This article delves into the methodology surrounding alterations in drop size and distribution along the column, elucidating their significance.In addition, methods based on artificial intelligence (AI) such as multilayer perceptron neural networks and gene expression programming were combined with a dimensional analysis approach to provide a new approach to estimating the Sauter mean drop diameter (d 32 ).In order to evaluate the proposed model, the data were obtained from the laboratory analysis.

Previous work
For all extraction columns, drop size is very important in hydrodynamic behavior and mass transfer performance.In other words, the size of the drops affects the presence and residence time of the dispersed phase, it affects the overflow and the interphase surface for mass transfer, and it affects the kinetics of mass transfer due to the mass transfer coefficients that are dependent on the drop diameter.In all extractors, the drop size distribution is not uniform along the length of the column, and the Sauter mean drop diameter is used as the average diameter of diameter size distribution index, which is determined by the relative rate of drop breakage and coalescence and is defined as follows: In pulsed columns, the dispersed phase is usually dispersed by a distributor, and usually, the drops have an average diameter between two and five millimeters (two for systems with very low interphase tension and five for systems with high interphase tension).
Kolmogorov 15,16 predicted the average drop size in pulsed sieve-plate extraction columns based on the theory of isotropic turbulence, which related the drop size to the input mechanical energy. (1) www.nature.com/scientificreports/Hinze 17 and Shinnar and Church 18 modified the theory of isentropic turbulence by means of an equation with the surface energy produced by the breaking of drops and the kinetic energy caused by the fluctuations of the same drop volume in the vortex on each side.
where the power loss per unit mass of liquids (ψ) can be calculated by the equation proposed by Jealous and Johnson 19 : For the orifice coefficient (C 0 ), a constant value of C 0 = 0.6 is recommended.The proposed equations for determining the average diameter of drops based on the Kolmogorov theory describe the effect of pulsation intensity on the size of the drops and predict the diameter of the drops in the range of the mixer-settler regime.Literature missing about the breaking of drops in pulsed sieve-plate extraction columns show that the breaking of drops when passing through the holes of the mesh tray is caused by their shear stress.Anyway, the theory of isotropic turbulence is not suitable for describing the process of drop breaking in mesh tray columns 20 .
Use On the one hand, new physical models have been proposed to predict drop size in pulsed sieve-plate extraction columns.Pietzsch and Pilhofer to model the stability of the drop when passing through the sieve-plate by equating the drag, buoyancy and inertia forces (which cause the drop to break) with the interphase tension force (which causes the drop to remain stable) gave a relation to have achieved 21 : A summary of other correlations to predict the Sauter mean drop diameter in pulsed columns is given in the Table 1. (2) Table 1.A summary of some predicted correlations of the Sauter mean drop diameter for pulsed columns.
Kerosene Toluene n-butyl acetate water Dilute sulfuric acid (4 g/L) water www.nature.com/scientificreports/Angelov et al. 14 obtained a relation based on computational fluid dynamics (CFD) for pulsed disc and doughnut columns when a drop breaks into two drops.
k is a constant depends on pulse intensity.
Tenova columns have been the subject of limited studies in the literature.Among the first studies reported by Li et al. 29 , the hydrodynamic performance of a Tenova pulsed column with a standard disc-doughnut pulsed column with a length of 200 cm and a diameter of 7.6 cm for copper extraction using LIX 84 compared.As a result, in the same operating conditions, less dispersed phase holdup was observed in Tenova column compared to disc and doughnut.At low pulsation intensity, the Sauter mean drop diameter was observed to be larger in the Tenova column compared to the disc-doughnut column.In another study by Li et al. 12 , the hydrodynamic performance of a Tenova column with a standard disc-doughnut pulsed column for the 336-Shellsol 2046 system was compared.Dispersed phase holdup and Sauter mean drop diameter under different pulsation intensities and different flow rates of phases.Their results showed that compared to a standard disc-doughnut pulsed column, the drops are larger and the dispersed phase holdup is lower.

Liquid-liquid system
The liquid-liquid system used in this study is the standard chemical system recommended by the European Federation of Chemical Engineering (E.F.C.E.).Toluene-water (high interfacial tension), isobutanol-water (low interfacial tension), and isobutyl acetate-water (medium interfacial tension) were used.
The physical properties of these systems, including the viscosity and density of each phase and surface tension, are shown in Table 2. Water saturated with solvent was used as the continuous phase and solvent saturated with water was used as the dispersed phase.Viscosity and surface tension were measured using an Ostwald viscometer (schott-Gerate GmbH, Germany), and a digital tensiometer (KRUSSN, Germany, model K10ST), respectively.

Equipment and experimental procedure
A schematic of the set-up used in this study is shown in Fig. 1.The active part was a tubular column with a diameter and length of 7.4 and 73 cm, which is made of Pyrex material.Its internal equipment consisted of 30 pairs of modified discs and doughnuts made of 0.2 cm thick stainless steel sheets which were placed alternately and at a distance of 1.2 cm from each other by means of Teflon spacers and were held by three support rods.The diameter of the discs, doughnuts and the holes of the doughnuts is 6.4, 7.4 and 3.4 cm, respectively and the free area fraction was 21.1%.Two precipitators were located at the top and bottom of the column with a diameter of 11.2 cm to completely separate the phases before exiting the column.
The drop sizes in this study were measured by photographing different parts of the column using a D3500 Nikon digital camera.The method of photographing the column was that after setting up the device in each test and setting the flow rate of the dispersed and continuous phase in the column, first the pulsation intensity was adjusted by the board controller and regulator, and the flow inlet valves were opened and then the pulse switch in the system turned on.After setting up the system, the aqueous and organic phases were allowed to flow in the column and about three volumes of the column, the flow of phases passed along the length of the column; In this case, the size of the drops remains constant and the drops reach their real size.At this time, when there is a pulse in the system, it was started to take pictures of different parts of the column.Photographs were taken from each stage separately along the length of the column and the photographs were used to calculate the diameter of the drops that were not spherical using Digimzer software.For elliptical drops both the minor and major axes, d 1 and d 2 were measured and the equivalent diameter, de, was calculated from Eq. 6, and an average of 200-300 drops were examined in each experiment and in total more than 10,000 drops have been analyzed.To examine the droplets, specific stages from the lower, middle, and upper segments of the column were individually selected under various test conditions.The software was employed to analyze the images captured from each stage.Calibration of the software relied on the thickness of discs and donuts as criteria, measuring all observed droplets in the images.Figure 2 provides actual images of the column as an example, showcasing analyzed images from stages 5, 7, 15, and 20, along with the respective count of analyzed droplets (5)  We c = 0.The physical characteristics of the systems used at a temperature of 22 °C are shown in Table 2.

Results and discussions
The effect of operational parameters on the Sauter mean drop diameter For all three systems at stage 2, the Sauter mean drop diameter decreases with increased pulsation intensity (see Fig. 3a).Similar results were obtained by the literature in this type of extraction column 29 .
It is expected that with the increase of the pulsation intensity, the Sauter mean drop diameter will reach its final limit, and the pulsation intensity at that point can be named the critical pulsation intensity, which can be the ultimate limit of breakage.
For all three systems, a decreasing trend is observed with increasing pulsation intensity in stages 2, 10 and 15 (Fig. 3a, b and c respectively).Comparing these three stages together shows that this reduction process starts from an average diameter of stage 10 less than stage 2 and in stage 15 an average diameter less than stage 10 and will reach the same final level at the point of critical impact intensity.
With the increase in the flow rate of the dispersed phase in stages 2, 10, and 15, as seen in Fig. 4, the Sauter mean drop diameter increases.It seems that with the increase in the flow rate of the dispersed phase, the average diameter produced in the distributor is larger until reaching the jet point in the distributor, so the increasing trend of the average diameter of the dispersed phase seems quite logical.It shows an increasing trend for three systems, from isobutanol to toluene, due to the increase in interphase tension.Similar results were obtained in the literature in disc and doughnut impact extraction columns 26 .
The decreasing trend of these changes from stage 2 to 15 for all three systems indicates that there is a breakage along the length of the column, so the decreasing trend of the Sauter mean drop diameter along the length of the column (from stage 2 to 15) seems quite logical and this trend has been seen in all three systems.
As can be seen in Fig. 5, the effect of continuous phase flow rate on the Sauter mean drop diameter is very small, but its increasing trend indicates the possibility of drop coalescence in different positions.
The decreasing trend of the Sauter mean drop diameter is in terms of the breaking of the drops along the length of the column, which indicates that the pulsation intensity was higher than the critical point ( Af cr ) and the breaking along the length of the column leads to a decrease in the Sauter mean drop diameter is from stage 2 to stage 15 for all three systems.

Drop size distribution
The drop size distribution for the toluene-water system and for stages 2-7 is shown in Fig. 6a.It is evident that the breakup process occurs during the column operation and leads to a tendency towards a normal distribution.
As can be seen, the number of drops with a larger size in stage 2 is much more than the drops with a smaller size, which indicates that the process of breaking is underway, and from stages 2 to 7, the distribution tends toward normal due to the breaking of the drops.It can be seen that the distribution will be normal from stage 7 onwards.This trend can be observed for the isobutyl acetate system as well (Fig. 6b), but as expected, the drop diameter is smaller than that of toluene.The trend mentioned for the isobutanol-water system can also be seen in Fig. 6c.It appears that along the length of the column, droplet breakage contributes to a normal distribution, predominantly observed in the middle stages.Following these middle stages, a consistent normal distribution is maintained due to a balance between breakage and coalescence, resulting in stability beyond this point.

The correlation of predicting the Sauter mean drop diameter
To achieve the design of such columns, we need correlation to predict the Sauter mean drop diameter.
As mentioned before, no correlation has been provided to predict this parameter specific to the Tenova column.In this study, a correlation based on operational parameters and physical characteristics of the system specific to this column is provided.In addition to predicting the Sauter mean drop diameter, changes in column length can be predicted.According to the background of the research and the identification of variables that affect Sauter mean drop diameter, using dimensional analysis, 9 products without independent dimensions were obtained to predict Sauter mean drop diameter.In order to create a linear equation between independent variables without dimension and d 32 , linear multivariate regression was used (see Figs. 7 and 8), and its equation is as follows: The constants of linear and non-linear regression are given in Table 3.
In this study, the gene expression programming (GEP) model was used to predict d 32 with high accuracy and also to validate the model with higher reliability.
The main advantage of this method is that, unlike other artificial intelligence methods, it provides the user with an equation.To model d 32 using this technique, the parameters of the model should be set first and then the model should be developed.
Here, 4 genes have been obtained, these genes must be connected to each other to create an equation.Therefore, the main equation used to predict d 32 is as follows: (8)   A comparison of the predicted and measured values using the GEP model is shown in Fig. 9.These results show that the value of d 32 is predicted with very high accuracy and the validation of the model is done with high confidence.
In order to compare the developed model with other artificial intelligence methods, a multilayer perceptron (MLP) artificial neural network was used as a powerful artificial intelligence method.This method has the ability (10) Two correlation for predicting the Sauter mean drop diameter reported in the literature 27,28 were evaluated for the experimental data of the present study.The average absolute relative deviation (AARD) was reported.
Figure 11 compares the Sauter mean drop diameter measured experimentally with the predicted value by the mentioned correlations.
In Table 4, the comparison of previously reported correlations for predicting the Sauter mean drop diameter and the present work.

Sensitivity analysis
Among the useful measures after modeling is determining the sensitivity of the target parameter to the input parameters.Typically, to determine the sensitivity value of the input parameters on the target, by removing one of the input parameters, the changes in the error value of the test data are checked.The severity of the increase in error indicates the greater impact of the omitted parameter on the results.In this study, after applying the sensitivity using the Cosine Amplitude Method (CAM), the intensity and impact of the input parameters on the output was calculated.This method was first developed by Yang and Zhang 30 and then by Jong and Lee 31 and used to determine the similarity of correlations between related parameters.In this method, an m-dimensional space where m is the number of input parameters is assumed: Each member of the input parameters such as X is connected to the objective function by a length vector: Therefore, each data point is a point in the m-dimensional space that requires m coordinate components for a complete description.The strength of the correlation (r) between both variables can be calculated as follows:  where x i and x j are independent and dependent variables, respectively.In this way, the impact of each of the input parameters is calculated.The greater the impact of the input parameter on the target (output), the r ij approaches 1, and on the other hand, if there is no impact of the input parameter on the target output, the value of r ij tends to zero.Normally, the value of r ij > 0.9 indicates a significant effect of the independent parameter on the output, and values < 0.8 indicate its weak effect on the output variable.Regarding the data used in this study, the sensitivity level of each of the input parameters is shown in Fig. 12.The results show that all independent parameters have an effect on d 32 changes.

Conclusions
The Tenova column, which has been introduced in recent years, has been built and tested on a pilot scale.Equations ( 8), ( 9), (10) were proposed for predicting of mean drop size along the column using the dimensional analysis approach.The proposed correlations are compared with the correlation presented in the literature.The results obtained from this work in the middle stages where the distribution of drops is completely normalized are in good agreement with the results of other researchers (error less than 29%) and an error of less than 12% has been obtained by editing the correlation presented with the GEP method.The correlation presented based on the stage number can be used in the length of the column, which is strongly recommended for the design of such columns instead of the average diameter.In addition to the average diameter and its changes along the length of the column, the distribution of drops and how it changes along the length of the column has been obtained experimentally for all the mentioned systems.The results are justifiable.For each distribution curve, at least 200-300 drops and in total, more than 10,000 drops have been analyzed in different operating conditions in this research work.

Forecasting 3 A
new type of Tenova pulsed extraction column was introduced in 2017.It is the newest generation of pulsed columns.Due to the internal equipment of this column and the lack of moving parts and the simplicity and speed of repairs and maintenance, it has been the focus of researchers in recent years.

Figure 1 .
Figure 1.Schematic of the set-up used in this study, Tenova pulsed column and internal equipment of the column.

Figure 2 .
Figure 2. Real images of the column to check the drop size along the column related to various stages, (a) 5th stage, (b) 7th stage, (c) 15th stage, (d) 20th stage.

Figure 3 .
Figure 3.The effect of pulsation intensity on the Sauter mean drop diameter for all three systems at (a) stage 2; (b) stage 10, and (c) stage 15.

Figure 4 .
Figure 4.The effect of dispersed phase flow rate on the Sauter mean drop diameter for all three systems at (a) stage 2; (b) stage 10, and (c) stage 15.

Figure 5 .
Figure 5.The effect of continuous phase flow rate on the Sauter mean drop diameter for all three systems at (a) stage 2; (b) stage 10, and (c) stage 15.

Figure 9 .Figure 10 .
Figure 9. Predicted and measured values of d 32 in the (a) training and (b) testing section of the GEP model.

Figure 11 .
Figure 11.Comparing the Sauter mean drop diameter experimentally with the predicted value; (a) Wang et al.(b) Sarkar et al.

Figure 12 .
Figure 12.Value of r ij factor in d 32 sensitivity analysis.

Table 2 .
26 = ρ c kd max σ Physical characteristics of the systems used at 22 °C.

Table 3 .
Constants of linear and non-linear regression.

Table 4 .
AARD for predicting the Sauter mean drop diameter based on the previous and the present work.